Correlations of Structure and Dynamics in an 
Aging Colloidal Glass 



Gianguido C. Cianci, Rachel E. Courtland and Eric R. Weeks* 

Department of Physics, Mail Stop 1131/002/1AB, Emory University, Atlanta, GA 

30322, U.S.A. 



Abstract 

We study concentrated colloidal suspensions, a model system which has a glass tran- 
sition. Samples in the glassy state show aging, in that the motion of the colloidal 
particles slows as the sample ages from an initial state. We study the relationship 
between the static structure and the slowing dynamics, using confocal microscopy 
to follow the three-dimensional motion of the particles. The structure is quantified 
by considering tetrahedra formed by quadruplets of neighboring particles. We find 
that while the sample clearly slows down during aging, the static properties as mea- 
sured by tetrahedral quantities do not vary. However, a weak correlation between 
tetrahedron shape and mobility is observed, suggesting that the structure facilitates 
the motion responsible for the sample aging. 
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1 Introduction 



When some liquids undergo a rapid temperature quench they can form glasses. 
This occurs at a glass transition temperature T g which often depends on 
the cooling rate. As the system is cooled, the approaching glass transition is 
marked by a dramatic increase in the macroscopic viscosity of the liquid, and 
a corresponding increase in the microscopic time scales for motion [1,2,3,4]. 
Both the viscosity and the microscopic relaxation time can change by many 
orders of magnitude as the temperature decreases by merely 10%. Once in the 
glass state, another phenomenon is noted, that of aging: the dependence of 
the properties of the system on the time elapsed since reaching T g . When such 



For correspondence: weeks@physics.emory.edu 



Preprint submitted to Elsevier Science 



6 February 2008 



behavior is observed the system is said to be out of equilibrium, a fact that 
could be anticipated by noting that the dependence of T g itself on cooling rate 
implies the glass transition is not an equilibrium phenomenon. Aging most 
prominently manifests itself in the dynamics: the microscopic relaxation time 
scale depends on the age of the sample. 

Attempts to explain these phenomena try to link the microscopic structure to 
the microscopic dynamics. For example, one might postulate that the increase 
in viscosity is caused by the growth of domains whose dynamics are correlated 
[5]. However, no experiment has seen a structural length scale characterizing 
such domains that grows or diverges at T g [1,6,7]. Likewise, aging might be due 
to some coarsening of structure; as domains of glassy structure grow, motion 
is slowed, which in turn slows the further growth of these domains. However, 
these domains have not been identified, and currently no structural features 
have been identified that explain aging dynamics [8,9]. 

Recently, some interesting developments in the study of aging non-equilibrium 
systems have been brought about by the adoption of dense colloidal suspen- 
sions as model systems for liquids, glasses and gels [7,8,10,11,12,13]. Colloidal 
suspensions consist of solid particles in a liquid, and the motion of the parti- 
cles is analogous to that of atoms or molecules in a more traditional material 
[10,14,15]. In these systems, the particle interactions can easily be tuned from 
repulsive to attractive. A common case is when particles interact simply as 
hard spheres with no interactions, attractive or repulsive, other than when 
they are in direct contact [10,16]. In all cases, a major control parameter is 
the packing fraction 0. For hard spheres this is the only control parameter 
and when <fi is raised above a value of g ~ 0.58 the system becomes glassy 
and the aging process begins. 

While structural changes in an aging system remain unclear, two experiments 
studying colloids have characterized the dynamics. Cipelletti and co-workers 
[12] studied aging in a colloidal gel using novel light scattering techniques and 
showed that the dynamics in such a non-equilibrium system present strik- 
ing temporal heterogeneities. Aging has also been studied in a colloidal glass 
using confocal microscopy [13]. In that study both temporal and spatial het- 
erogeneities were seen. However, despite the ease with which a colloidal glass 
can be formed and observed a detailed understanding of the structural changes 
that accompany aging and the slowing of dynamics has not yet been seen. 

In this paper we study the structure of an aging colloidal glass by considering 
how colloids pack together. Entropy can be maximized by optimizing packing 
in a dense suspension. Consider the intriguing case of the crystallization of 
hard spheres. When spheres arrange into a crystalline lattice, they lose con- 
figurational entropy. However, they each have more local room to move close 
to their lattice site, and thus the vibrational entropy is larger. This increase 
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in vibrational entropy outweighs the loss of configurational entropy due to 
crystallization [17]. In practice, this argument holds true for systems with vol- 
ume fractions above 0f reez e = 0.494, the point at which the system begins to 
nucleate crystals; below 0f ree ze, the configurational entropy dominates and the 
system prefers an amorphous, liquid configuration [10]. 

For glasses, we consider a different sort of packing. The most efficient way 
to pack four spheres of diameter d in three dimensions is to place them at 
the vertices of a regular tetrahedron with edge length d. In this configuration, 
the effective volume fraction for the four spheres reaches a surprising 0.78. In 
other words, for a given volume fraction 0, locally four particles can maximize 
their entropy by arranging into a regular tetrahedron consistent with the global 
volume fraction <p, thus giving them additional room to move. However, regular 
tetrahedra do not tile 3D space and therefore the most efficient macroscopic 
packing is that of a hexagonally packed crystal at 0hc P ~ 0.74. Thus in a 
glass there is a frustration between the drive to locally pack in tetrahedra to 
maximize the local volume available to vibrations, and the inability to tile 3D 
space with such structures. This has been suggested as a possible origin for 
the glass transition in simple liquids [18,19,20,21]. 

We take advantage of the insight afforded to us by fast laser scanning confocal 
microscopy [7,14,15,22] and study an aging colloidal glass in terms of tetrahe- 
dral packing. We focus on geometrical properties of the tetrahedra formed by 
the colloids and look for correlations between these static quantities and the 
conspicuous slowing of the dynamics as measure by the average tetrahedral 
mobility. We find that while the distribution of these static quantities does 
not age, they correlate weakly with mobility, suggesting that the structure 
facilitates the aging process. 



2 Experimental Methods 



We suspend poly(methyl methacrylate) (PMMA) colloids of diameter d = 
2.36yum in a mixture of 15% decalin and 85% cyclohexylbromide by weight. 
The mixture closely matches the density and refractive index of the parti- 
cles, thus greatly reducing sedimentation and scattering effects. The size poly- 
dispersity of the colloids (~ 5%) prevents crystallization. The particles are 
sterically stabilized against van der Waals attractions by a thin layer of poly- 
12-hydroxystearic acid [23]. We dye the colloids with rhodamine 6G [22]. The 
particles also carry a slight charge due to the dye. In this paper, we measure 
all lengths in terms of the diameter d and all times in terms of Tdiff, the time 
a particle would take to diffuse its own diameter in the dilute limit. Given the 
solvent viscosity rj = 2.25 mPa-s at T = 295 K, this time is = 11.4 s where 

£) = k B T 
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We acquire three dimensional images by fast laser scanning confocal mi- 
croscopy at a rate of 1 every 26 s. The observation volume measures 26d 
x 25d x 4.2d. At these high densities (0 > 0.58) the colloids move slowly and 
can easily be tracked using established analysis techniques [15,22,24]. The 3D 
positions of ~2500 particles are measured virtually instantaneously with an 
accuracy of 0.013<i in the x-y plane and 0.021<i along the optical axis. We 
acquire data at least 25<i away from the closest wall, to avoid boundary effects 
[25,26]. We do not observe any crystals in the bulk even after several weeks. 

The phase behavior of this quasi-hard sphere system is controlled by varying 
the packing fraction <fi. The system undergoes a glass transition when <fi > 
(f) g ~ 0.58 in agreement with what is seen in hard sphere systems [10,15]. Here 
we present data from a sample at ~ 0.62 though we see qualitatively similar 
results for all > <p g . 

Proper sample initialization is paramount when studying aging and is ensured 
here by a vigorous, macroscopic stirring. This shear melting effectively reju- 
venates the glass and yields reproducible dynamics that depend exclusively 
on t w , the time elapsed since initialization. Data acquisition starts immedi- 
ately after rejuvenation. Transient macroscopic flows are observable for the 
first 25 min^ 140r d iff and we set t w = 0, or age zero, when they subside. The 
results below are insensitive to small variations in this choice. 



3 Results 

We observe our sample for ~ 700r d ig- without disturbing it. We then split 
the data in three time windows as follows: [0 — 100r diff ], [100 — 300r diff ] and 
[300 — 700r d ifr]. This corresponds to doing three experiments with samples 
aged t w = 0, 100 and 300r diff respectively. The dynamics slow as the sample 
ages, as shown in Fig. 1, where we plot the mean square displacement for the 
three data portions averaged over all particles and over all initial times within 
a given time window. At short and medium times < 10), particle motions 
are subdiffusive as indicated by a slope less than unity on the log-log plot. At 
longer times the slope tends to one; the time scale for this upturn changes 
dramatically for different values of age t w clearly indicating that the sample is 
out of equilibrium. It is this slowing down of dynamics that we wish to analyze 
in terms of tetrahedral structure. 

We start our structural analysis by calculating the pair correlation function 
g(r) and plot the result in Fig. 2. This function does not vary with age and thus 
is calculated by averaging over all times. The first peak of g(r) is at r = 1.04c? 
which deviates somewhat from the ideal hard-sphere position (r = d). This 
can be explained by the slight charging mentioned above and perhaps also by 
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Fig. 1. Aging mean squared displacement for a colloidal glass at <j> ~ 0.62. The three 
curves represent three different ages of the sample. A : t w = Ordifr, x : i w = lOOrdiff 
and O : *w = 300rdiff. The dashed line has a slope of 1 and represents diffusive 
behavior, not seen in this glassy sample. 

the uncertainty in the value of the particle diameter which we deem to be at 
most 2%. Figure 2 also shows the characteristic double second peak found in 
many glassy systems. 

In order to study the tetrahedral packing in our sample we begin by labeling 
as nearest neighbors every pair of colloids whose separation is within the first 
peak of g(r), namely 0.7 Ad > r > 1.38c?, as is shown by the shaded area 
in Fig. 2. The lower limit is chosen to eliminate artificially close pairs which 
arise from the occasional error in particle identification, while the upper limit 
corresponds to the first minimum of g(r). Note that a completely coplanar 
arrangement of four spheres in a square is excluded as a tetrahedron, as the 
diagonal would have length y/2d which is excluded by our upper limit. The 
results presented here are insensitive to small variations in these parameters 
and match those obtained using Delaunay triangulation as a nearest neighbor 
finding algorithm. 
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Fig. 2. Pair correlation function g(r). The shaded area indicates the range of inter- 
particle distances used to define nearest neighbors. 
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A tetrahedron is then defined as a quadruplet of mutually nearest neighbor 
colloids. To characterize each tetrahedron, we calculate several geometrical 
characteristics. The first is "looseness" b, defined as the average of the lengths 
of the 6 edges, or "bond lengths", 6;. An "irregularity" a b is defined as the 
standard deviation of the b\. The looseness and irregularity appear to be the 
most important geometric parameters to characterize a tetrahedron shape, as 
will be discussed below. The nondimensional volume V/ d 3 and nondimensional 
surface area A/d 2 are also measured. To quantify an effective aspect ratio of 
each tetrahedron, we calculate the height of the tetrahedron as measured from 
each of the four faces, and consider the largest height H and shortest height 
h. We form aspect ratios from these two heights by dividing by the areas 
(Ah and Ah respectively) of the tetrahedron face they are perpendicular to. 
Conceptually, large values of H 2 /A H correspond to thin pointy tetrahedra, 
and small values of h 2 /A h correspond to flat pancake- like tetrahedra. We thus 
term these two quantities "sharpness" and "flatness" respectively. 

In addition to these structural characteristics, we also consider the dynamics 
by the tetrahedral mobility /i, which is calculated by averaging the distances 
moved by the four colloids over a time At = 50r diff : 

M*) = 2E|Ar1(t,Af)| (1) 



The results that follow do not depend sensitively to the choice of At. 

To assess the value of these structural and dynamical characteristics, we cal- 
culate the correlation coefficients between \i and the other tetrahedral char- 
acteristics and show them in Table 1. This is done in the standard way of 
defining correlation coefficients, 

_ l ^ (gi - v) (ft - q) (0 , 
Jy i=1 o v (J q 



where p and q are any two variables with averages p and q, and standard 
deviations a p and a q . In our case the sum runs over all tetrahedra and all 
times. In Table 1, a value of one would signify perfect correlation, a value of -1 
would represent perfect anti-correlation while a value of zero would indicate 
completely uncorrelated data. 

Given that we are trying to understand the slowing of the dynamics seen in 
Fig. 1, we focus on the correlation between mobility fi/d and the structural 
characteristics. While all the coefficients are quite small, those that relate 
mobility to looseness b/d and irregularity cr b /d are relatively big. Mobility is 
also noticeably anticorrelated with the flatness h 2 /A h . Some insight into these 
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Table 1 

Correlation matrix for some geometrical characteristics of tetrahedra and mobil- 
ity. The matrix is symmetric with respect to the diagonal so the lower half is not 
repeated, b is the average length of the tetrahedra edges ("bonds") and a b is the 
standard deviation of these lengths. See text for details of the other characteristics. 



correlations is gained by considering the correlations between the different 
structural characteristics, as shown by the other entries in Table 1. The flat- 
ness h 2 1 'Ah is strongly anticorrelated with irregularity a b , and given the more 
intuitive nature of a b and its simpler mathematical definition, in what follows 
we focus on a b rather than h 2 /A h . The volume and area parameters, V/d 3 and 
A/d 2 , are strongly correlated with the looseness, which is sensible given that 
they all measure the size of a tetrahedron. 

We therefore choose to study the looseness and irregularity as being both 
relatively well-correlated with mobility, both easily defined in terms of the six 
tetrahedron edge lengths, and weakly correlated with each other (as seen in 
Table 1). The last point suggests that they capture two distinct properties of 
tetrahedron structure which are both important for mobility. 

Figure 3 shows the distribution of tetrahedra in the b, ovplane. The closed 
curves represent the levels of abundance of tetrahedra with a given value of 
b and a b with respect to the abundance of the most probable tetrahedron at 
b ps lAld and a b ps 0A2d. Somewhat surprisingly, the distribution does not 
age [27] and we therefore take an average over all times. Figure 3 shows a 
broad variability of both looseness and irregularity. However, these curves do 
outline a major axis along which many of the tetrahedra lie. This axis suggests 
that the looser the tetrahedron the more irregular it is bound to be. This is 
reflected in the correlation coefficient of b and a b in Table. 1, although its 
relatively small value (0.37) highlights the breadth of the overall distribution. 

We present a qualitative picture of the relationship between the above two 
static geometrical quantities and the dynamic quantity of mobility in Fig. 4 
which shows the average value of mobility as a shade of grey. The more mo- 
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Fig. 3. Contour plot showing the abundance of tetrahedra with a given looseness and 
irregularity. The iso-curves are labeled relative to the peak tetrahedral abundance 
at {b/d= 1.11; a b /d = 0.12}. 
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Fig. 4. Plot of tetrahedral mobility (^(At = 50^^)), averaged over all ages, versus 
looseness b/d and irregularity a b /d. The darker the color the more mobile the tetra- 
hedron. The contour lines are the same as in Fig. 3 and represent the abundance of 
tetrahedra with a given value of b/d and a b /d. 

bile combinations of b and a b are darker and the less mobile combinations are 
lighter. This rendering gives a clear qualitative view of the correlations be- 
tween these quantities. Specifically, we note that mobility increases with both 
looseness and irregularity of the tetrahedron. This makes intuitive sense: a 
larger value of b suggests a smaller local volume fraction and thus more room 
to move, and a larger value of a b likewise suggests a poorly-packed structure 
which may have more room to move. This also agrees with previous results 
seen in supercooled colloidal fluids [28,29,30]. Overplotted on the intensity plot 
of Fig. 4 are the same abundance contours as seen in Fig. 3, showing us that 
the most probable tetrahedra are a medium shade of grey - they are neither 
the fastest or slowest tetrahedra. 

We thus have two results, an overall slowing of dynamics seen in Fig. 1, and 
a relationship between structure and dynamics seen in Fig. 4. This suggests a 
possible hypothesis for aging, that the slowing of the dynamics is an accumu- 
lation of structure corresponding to slower dynamics: a buildup of tetrahedra 
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Fig. 5. Average tetrahedral mobility as a function of looseness b/d. The three curves 
represent three different ages of the sample. A : t w = Cbdifr, x : i w = lOOrdifr and 
O : = 300rdifr. The dashed curve represents the distribution P(b/d) and is shown 
here to highlight the lack of statistics at low (b/d < 0.98) and high (b/d > 1.2) values 
of b/d. 

with small values for looseness and irregularity. As mentioned earlier, though, 
the overall distribution of tetrahedra structural properties (Fig. 3) does not 
depend on the age of the sample [27]. To reconcile this, we consider in more 
detail the connection between structure, dynamics, and the age of the sample. 

We show the influence of looseness on the tetrahedron mobility in Fig. 5 where 
we plot the average mobility of tetrahedra as a function of looseness. We do 
so averaging over the three sample ages separately. If we consider each curve 
separately, we note a reproducible trend: the least mobile tetrahedra are those 
with b ~ d indicating that those tetrahedra are very tightly packed. At very 
low values of b mobility increases somewhat, although as mentioned above, 
tetrahedra with b/d < 1 may be erroneous. Furthermore, note that there 
are extremely few tetrahedra with b/d < 1 as shown by the dashed curve 
representing P(b/d), the probability of finding a tetrahedron with a given 
looseness averaged over all ages. We therefore cannot put too much weight on 
the values of /x for b < d. (The same can be said for tetrahedra with b > 1.17c?.) 
In the intermediate range there is a clear trend that looser tetrahedra are 
more mobile. Thus the structure in some way facilitates the aging, in that 
looser regions are more free to rearrange. However, it is also important to 
note that each symbol is an average over a broad distribution of mobilities 
associated with the given value of b/d. In particular the standard deviation 
of the distribution is almost comparable with the average value. This simply 
means that the correlation between b and /i is a weak, average effect and not, 
for example, a usefully predictive relationship [30]. 

As expected with any plot involving the mobility of this system, aging is clearly 
visible in Fig. 5 as the three curves are shifted down to lower mobilities as t w 
increases. The overall shape of the curves, however, does not depend on the 
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age of the system. In other words, we are not witnessing a relative shift in the 
mobility of tetrahedra with varying looseness but merely an overall slowing 
down of all tetrahedra. 
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Fig. 6. Average tetrahedral mobility as a function of irregularity a^/d. The three 
curves represent three different ages of the sample. A : t w = Ordifr, x : t w = lOOTdiff 
and O : = 300rdiff. The dashed curve represents the distribution P{ab/d) and 
is shown here to highlight the lack of statistics at low (a^/d < 0.01) and high 
((Jb/d > 0.2) values of Objd. 



A similar analysis on the relation between tetrahedral mobility and irregularity 
is shown on Fig. 6. As a reference we plot P(<Jb/d), the probability of finding 
a tetrehedron with a given irregularity. This distribution does not age and 
is therefore averaged over all times. Again, we look at the dependence of /i 
on Ob for the three ages in our experiments and find that there is a positive 
correlation, as previously indicated in Table 1. Just as in the case of b, this 
positive correlation is an average effect and again, the distribution that leads 
to each of the symbols on the figure is quite broad. Nevertheless, a reproducible 
difference of ~ 10% in the mobility differentiates very regular tetrahedra from 
very irregular ones. Again, just as above, aging is evident in the data and 
again it has no strong effect on the shape of the curves but rather uniformly 
slows down tetrahedra with all values of irregularity. 



4 Conclusion 



We observe colloidal glasses and find clear signs of aging in the mean squared 
displacement of the particles (Fig. 1). We analyze the static structure of the 
aging sample in terms of tetrahedral packing of colloidal particles. We find a 
broad distribution of tetrahedra as measured by the distributions of tetrahe- 
dral "looseness" and "irregularity" , corresponding to the tetrahedron's mean 
edge length and the standard deviation of edge lengths, respectively (Fig. 3). 
These two quantities are slightly correlated; on average, the looser a tetrahe- 
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dron is the more irregular it will be. More importantly, we find that tetrahedral 
shape and mobility are somewhat correlated: the looser and the more irregular 
a tetrahedron is the higher its mobility (Fig. 4). This suggests that aging might 
be due to an increase in tight, regular tetrahedral structure, but surprisingly 
the distribution of geometrical quantities is age-independent. Instead we find 
that aging indiscriminately affects tetrahedra with all values of looseness and 
irregularity by uniformly decreasing their mobility. 

In conclusion, we find that static structure as measured by tetrahedral quan- 
tities does not indicate the age of a glass. None of the distributions of the 
geometrical quantities considered in Table 1 show any aging. However, at 
any instant in time the age of our sample must somehow be encoded in the 
positions of the colloids and, for example, analyzing the spatial correlations 
between tetrahedra, while beyond the scope of this paper, might provide more 
insight into this matter. Finally, since aging ought to result in subtle config- 
uration changes, and since the looser and more irregular tetrahedra allow for 
the most motion to happen, we can infer that the local structure does indeed 
facilitate the aging process. While it is worth noting that it is not established 
that the most mobile particles are the most important ones for aging, the 
connection between structure and mobility holds true for less mobile particles 
as well. Table 1 suggests that in this respect, tetrahedral irregularity is the 
most significant quantity whose positive correlation with mobility lends sup- 
port to our original motivating idea that perfect tetrahedra are an important 
structural element in a glass. 
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